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We have studied the phenomenon of weak localization in multi-quantum-well (MQW) structures 
in the regime of weak tunneling, when superlattice minibands are not formed. We have calculated 
the effect of weak localization on conductivity, which in this situation is described by a system of 
coupled Dyson equations. The tunneling across the MQW structure is found, in general, to suppress 
the weak localization effect on conductivity. 
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The weak localization, which results from interference 
of two electron waves propagating along a closed path 
in opposite directions and manifests itself, among other 
effects, in the phenomenon of the negative magnetoresis- 
tance, has been studied extensively in metals and semi- 
conductors, aSpWjSjll as thin films and quantum wells of 
these materialstJij. Recently, a lot of attention has been 
devoted to studies of the weak localization in more com- 
plex systems, such as suaeslattices, multiple quantum 
wells, quantum wires, etcclQ 

The physics of the weak localization in the su- 
perlattices is of particular interest, and has -been 
fairly extensively investigated, both theoretically!] and 
experimentallyfl. The theories vary in their approach, 
however, they all start from (implicit or explicit) assump- 
tion that the electron spectrum consists of well-defined 
minibands. This supposition is in fact equivalent to as- 
suming that the tunneling time between two wells t±2 is 
much shorter than the momentum lifetime for in-plane 
motion, t\\. This is a necessary condition for miniband 
formation, and for the structure to be treated as a super- 
lattice. 

The opposite case, when t\\ -C t\% and an electron 
diffuses in an single quantum well for a long time before 
tunneling into a neighboring well, has never been stud- 
ied, to the best of our knowledge. However, in practice 
both cases can be easily implemented, the first one - in a 
superlattice, the second one - in weakly coupled multiple 
quantum wells, fn this latter case, there are no mini- 
bands, and no coherent motion across the planes of the 
MQW structure. The tunneling is sequential, i.e. there 
is no tunneling through many barriers at once, and tun- 
neling into the next well does not depend on the history 
of electron motion. 

We will be interested in the effects of such weak tunnel- 
ing on the weak localization in lateral transport. Conse- 
quently, we choose the structure and make assumptions 
which help us to investigate this effect without additional 
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FIG. 1. Cooperon for multi-quantum-well structure, 
indices i, j, k, and I denoting quantum well numbers. 
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complications. For example, we take the simple model 
for electron spin relaxation, which is described by the 
single spin-relaxation time rgo = T Smx = i~ Syy = 2t Si u. 
From Refs. [ ||§ one can see that in A 3 B 5 quantum 
wells this approximation works quite well for large elec- 
tron densities. We also assume that the temperature is 
low, r < 6p, and that Ti/thEy <C 1 -rfhe usual condi- 
tions in the theory of weak localizationai (here e-p is the 
electron Fermi energy). 

The weak localization contribution to the conductiv- 
ity of a sipjjjc quantum well is given by the well-known 
expressio: 
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where 

C(q) Q /?/3a is the Cooperon, a and j3 are the spin indices 
which we will omit in all subsequent expressions, g^ ax = 
(-Dii) -1 , D = v 2 ti/2 is the diffusion coefficient, t\ is the 
transport time for in-plane motion in a well (differs from 
the lifetime Tn for long-range scattering), and vq = m/2ir 
is the density of states at the Fermi level. 

In the case of a MQW structure, the Cooperon is a 
4-dimensional matrix indexed by the well numbers, Cy 



(see Fig. y 
N wells is 



The total conductivity of the structure with 
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where C„„ = C"™. For large number of wells N the cal- 
culations can be greatly simplified if we impose periodical 
boundary conditions, so the TV— th well is connected with 
the 1st. We will also assume all wells to be identical. In 
this case 



Act 



TVActi 
e 2 D 



(3) 



«/3 n 1 



(2tt)2 



The Cooperon can be found from the Dyson equation, 
or, in our case, system of Dyson equations, since the com- 
ponent C nm is connected with C nm ±i. These equations 
are represented by diagrams shown in Fig. |[ Note that 
each "component" is actually a 4-dimensional 2x2x2x2 
matrix with spin indices, which can alternatively be rep- 
resented as a 4 x 4 matrix indexed by pairs a/3. These 
equations contain only components Ci„ = C}™; similar 
equations exist for components €2™ and so on. There are, 
of course, other components of the Cooperon CfJ , how- 
ever, none of them is relevant for the weak localization 
conductivity correction: some turn to after averaging 
over impurities (like C 2 \ which contains (V11V21) = 0, 
where Vy is the matrix element for scattering with ini- 
tial state in well i and final state in well j); other, while 
non-zero, do not contribute to to the conductivity, (for 
example, the contribution of C 22 is proportional to the 
overlap of the wave functions in different wells, because 
lines 1 and 2 would have to meet in the vortex, and 
we assume that this overlap is extremely small). These 
ignored diagram also does not contribute in the Dyson 
equations, unless the scattering in the two wells is corre- 
lated, i.e. they do not couple with Cn and C12 as long as 
(V11V22) = 0- Another component, C^, in principle can 
give contribution to the conductivity (just like its spin 
analog, C^, does give such contribution). However, in 
order for this component to be non-zero, there has to 
be at least one scattering act which occurs on the same 
scattering center but in different wells. We assume that 
the scattering centers have short enough range so there is 
no correlation in the scattering potentials in neighboring 
wells. Hence, there exists no such center which can scat- 
ter electrons in both wells, and the Cooperon component 
C 2 2 vanishes. 

Written as a system of integral equations, Fig. |^ be- 
comes 

Cn(q) = \Vn\ 2 
+ J d 2 9l |Vii|M(w,gi + q3Gr(w,-gi)C u (q) 

+ / d 2 g 2 |^ 12 | 2 G+( W ,g2 + q)G 2 -( W ,- g 2)C 2 i(q) 



+ / ' d 2 g N \V 1N \ 2 G%(uj,g N + q)G N (u,-g N )C N1 (q), 

C 2 i(q) = \V 2 i\ 2 (4) 
+ / d 2 9l |F 21 | 2 G+(o;,g 1 +q)Gr(w,-g 1 )Cii(q) 

+ J d 2 g 2 |^2| 2 G+(^,g 2 + q)G-(^,-g2)C 21 (q) 

+ / d 2 g 3 |y 23 | 2 G+( W ,g3 + q)G3( W ,-g 3 )C3i(q), 

and so on for all the components C n i- Note that the 
tunneling is only possible between the neighboring wells, 
so only terms with V n . n ±i are present. Here G ± (w,k) 
are the Green's functions, 



G ± (w,k) = 



1 
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E(k) = k 2 /2m and r v is the phase relaxation time. After 
the integration by E(g) in the right-hand side of Eq. (Q) 
the result is expanded up to second order terms in se- 
ries in small parameters th/t v and vqrn (v = dE /9k). 



Also, the Cooperon Cf 7 -(q) is expanded in harmonics: 



Cg(q) cosn0 q , 
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where tan</> q = q v /q x - These transformations follow 
closely Refs. J_0-f| and are described in details, for exam- 
ple, in Ref. [ Only the 0-th harmonic of the Cooperon, 
C^J = C%j , gives non-negligible contribution to the 
conductivity. For the components C n \ we arrive to the 
following system of linear equations: 

Gn = \Vn | 2 + — (1-£iito) Gn 

Til 

+ JH (1 _ L12T-0) (G12 + c N1 ) , 
T12 

G21 = |y 2 i| 2 + — (l-L n T )C 21 

Til 

+ — (1 - Li 2 r ) (Gn + G31) , (7) 
T12 

G31 = (1 — LuTq) G31 

Til 

+ ^-(l-L 12 T ) (G21+G41),... 

T12 

Here 
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T11 is the momentum lifetime in a well, T12 is a tunneling 
time between two neighboring wells: 



= 2^ \Vi n \ 2 , 
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and 
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is a 4 x 4 matrix, indexed, like the Cooperon itself, by 
pairs of spin indices a/3. In the basis of eigenfunctions of 
total momentum of the pair of electrons, singlet <po and 
triplet <j>™ with m = — 1, 0, 1 this matrix becomes 

L = Dq 2 + 1, L? = Dq 2 + ^ + i±^°, (9) 

1~4> T<j, T SO 

where S m o is the J-symbol. 

The system (Jjj) can be further simplified if we exploit 
the condition that ria 3> Tii ~ th, and keep only terms 
of the first order in tu/tis (however, the product Lt±2 
can be large or small, and we make no assumptions about 
its magnitude here). We can also simplify the expressions 
by introducing 

S(q) = 2™ T 1 2 1 C(g) (io) 
Then the system of linear equations for Si n becomes 



where 



AlSn + ^2 (521 + Sni) = Til, 

Ai5 3 i + A 2 (Su + Su) = Tia « 0, 

^l5 3 l+^2(521+5 4 l) = 0, ... 
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The system of linear equations ( |11[ ) allows us to find 
Sii(g) (which is a 4 x 4 matrix in the basis of singlet 
and triplet eigenfunctions <f>o, 4>T)- The weak localiza- 
tion correction to the conductivity (0) can be expressed 



through-ntbe-ieigenvalues of the inverse matrix, S 1 i(q), as 
followsBtia: 
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Here eigenvalue £o corresponds to the eigenfunction <pQ, 
and E™ - to 0™. Since the matrix L, and, hence, Su, 
are diagonal in this basis, the eigenvalues Eq and E™ 
of Sii (?) are just inverse of its corresponding diagonal 
elements, (S l n) and (Su)™, respectively, and 
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Here we have omitted indices 11 of S to simplify the 
notation. 

In a magnetic field B perpendicular to the planes of 
MQW the wave vector q becomes an operator with the 
commutation relations 
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where q± = q x ± iq y . This allows us to introduce cre- 
ation and annihilation operators a' and a, respectively, 
for which [aa^] = 1: 

D 1 ' 2 q + = t^V D x l 2 q^ = 5 l ' 2 a\ Dq 2 = 5{aa?}. 

(16) 

The non-zero matrix elements of these operators are 



(n —lan) = (n\ a T \n — 1) = vn, 



(n| {aa^} \n) = n + — . 
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The integration over o in a magnetic field becomes sum- 
mation over n, and the weak localization correction to 
the conductivity can be written as 
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where n max — 1/5th. It is convenient to extend the 
summation in Eq. (|l8| ) to n — > oo. One can do it using 
the following transformation, which exploits the condi- 
tion n max ^ 1 (this is the necessary condition for the dif- 
fusion approximation to work in the first place). We will 
see below that at large n the asymptotic values of So n and 
S"n are 1/Sn. Therefore, the expression under the sum in 
Eq. ( |lg| ) falls off as 2/5n and the sum diverges at large n. 
To remove this divergence we add and subtract 2/n+l to 

each term in Eq. ©: *E^" (-<%„ + Sft) = 

The first sum can be extended to n = oo because 
the expression under the sum falls of as 1/n 2 at large 
n. The second sum, again for n max 3> 1 can be ap- 
proximated by lnn max . Since the quantity of practical 
interest is not Act itself but the magnetoconductivity 
Act(£?) — Act(O), which is measured experimentally, we 
can replace lnn max = lnl/(5rn by Iii1/5t v and ignore 
the B- independent term In t v /th. After these transfor- 
mations, we arrive to the expression for the magnetocon- 
ductivity which does not contain t%i at all: 



Aa(B) - Act(O) = 

e 2 

- 21n(^)] 



(19) 
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The linear system ( |IT| ) with coefficients Eq. ( |i"2| ) re- 
mains the same in a magnetic field, with L now being a 
function of n instead of q: 
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The solution of the system (11) can be written in the 



following form, for each n and each spin index Z, m: 

1 



S = Sii = —F N (x) 



(21) 



where N is the number of wells in the MQW structure 
and x = Lt\2- The function Fn{x) can be found by 
solving the system (^lj) for any particular N. One can 
check by direct substitution that the following general 
expression is valid for any N: 



F N (x) = l(-x) N /\F 1 (^- 
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where 2-Fi(a, b; c; z) is the Hypergeometric function. For 
practical purposes it is more convenient to represent 
Fjq{x) by a ratio of two polynomials: 
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It is interesting to note that this expression has a very 
compact limit at N — > 00. By expanding F^(l/y) for 
very large N in series around y — one can see that 
the first TV coefficients of the series do not change with 
increasing N . They are given by the following expression: 
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To obtain the limit F oc (x) we extend the summation to 
infinity, and the sum can be written in a closed form: 
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In the above derivation of F^ (x) we have assumed the 
periodical boundary conditions, i.e. that electron can 
tunnel between the 1st and iVth wells. While greatly 
simplifies the calculations, this condition becomes inaccu- 
rate for small N, and a more accurate solution is needed. 
If there is no periodical boundary conditions, each well 
gives a different contribution into the total conductiv- 
ity (see Eq. ||)). The components of the cooperon C nn 
can be found from N systems of linear equations, system 
number m is written for coefficients C\ m , C^m, ■■■ , CWm, 
and is similar to Eq. ^ but does not contain the tunneling 
between 1st and JVth wells. From these systems we can 
find the components of the Cooperon C nn , or S nn , the 
sum of which enters the expression for conductivity. We 
can define an "average" to be used instead of Sn in 
the above expressions for conductivity, and, similarly to 
Eq. (^i"|), introduce function F^(x): 
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Figure [| shows the function F^ (x) for different num- 
ber of wells N, from 2 to infinity. For large N, the dif- 
ference between F^(x) and -Fjv(x) is negligible, and the 



4 



0.8 - 



0.6 



J3- 



0.4 



0.2 



FIG. 




Function 



determines Cooperon 



components entering the expressions for conductivity (see 
Eq. @) for N = 2, 3, 4, 6, and N = oo (solid curves, 
from highest to lowest). All functions shown by solid curves 
are obtained without using periodical boundary conditions. 
Using these boundary conditions for N = 4 results in Fi(x) 
shown by short-dash line. Note that the periodical boundary 
conditions become much more accurate if the x axis is scaled 
by (N — l)/N (for the case of N = 4, shown by long-dashes 
line, the abscissae axis was scaled by 3/4. 



two coincide exactly for N = oo. However, for small N, 
the periodical boundary conditions introduce a notice- 
able error; for example, in Fig. ^ we compare the func- 
tions F£(x) and F^(x) (short-dashed curve). It is inter- 
esting to note that replacing ri2 by an "effective" tunnel- 
ing time r{ 2 = N/(N-l)-rt2 (or xbyx' = N/(N-l)-x) 
makes the function Fn(x') almost coincide with F^(x). 
The reason for this behaviors is intuitively clear: by ap- 
plying periodic boundary conditions we have added an 
"extra" tunneling link to the structure, in addition to 
the N — 1 "real" links. Scaling T12 makes the "cumu- 
lative" tunneling rate for the structure with periodical 
boundary conditions equal to that of the real structure. 
For N = 2 this procedure, obviously, gives exactly cor- 
rect function, since the periodic boundary conditions at 
N = 2 are equivalent to just doubling the tunneling rate. 
In Fig. U we show by the long-dashed curve the scaled 
function i<4(4/3x); one can see that it looks very similar 



to Fl(x). 

We now have 
Eqs. ([19, 
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all the necessary expressions, namely, 
23) , to calculate the magnetoconductiv- 
ity due to the weak localization in an MQW structure. 
First, following— L[Il],|l2j we introduce the characteristical 
magnetic fieldsllij 
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and dimensionless conductivity -per one quantum well: 
2nh 



5a(B) 
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[Acr(B) - Act(O)] 
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We begin by analyzing two limiting cases of weak and 
strong tunneling, or x S> 1 and 1 « 1, respectively. In 
the first case we can immediately see from Eq. ( p4j ) that 
.Fjv(oo) = 1, so S = 1/L (note that at large n we also 
have x 3> 1 and, therefore, S = 1/L ~ 1/Sn). Sub- 
stituting expressions for L from Eq. (EO) we obtain the 
standard result for single quantum weho; 
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where $ is the digamma-function. This is not surprising 
since x — > 00 corresponds to T12 — > 00 in which case 
the parallel wells become completely independent. In 
the opposite case, x — > 0, we have from Eq. ( p3| ) that 
Fn(0) = 1/iV, so NS = 1/L and we again obtain the 
classical result, only now for the entire MQW structure 
since the tunneling couples the wells so that they all act 
as one well. 

We now present the results of our calculations of the 
weak localization correction to the conductivity in mul- 
tiple quantum wells for various parameters: number of 
wells N, spin relaxation time rso, and tunneling time 
Ti2- We begin by analyzing the dependence of the di- 
mensional magnetoconductivity da(B) on the number 
of wells N, which is shown in Fig. || for Hso/ Hp = 2 
this value is typical for A3B5 quantum wells, see Ref. [ 
Hu/Hp = 4, and Af from 2 to infinity. This figure 
also shows the effect of periodical boundary conditions at 
small A^: the magnetoconductivity calculated with such 
boundary conditions for A^ = 4 is shown by a dashed 
curve. One can see that increasing the number of wells 
makes the minimum on the .magnetoconductivity curve 
less pronounced. It is knowntlfEj that increasing rso has 
similar effect on the magnetoconductivity, the minimum 
becomes less pronounced and eventually disappears as 
rso increases; we illustrate this effect in the insert of 
Fig. I where Scr(B) is shown for N = 2, H 12 /Hp = 4, 
and Hso /Hp — 1/2 (top curve), Hso /Hp = 2 (middle 
curve), and Hso/ Hp = 3 (the lowest curve). One can 
say that coupling the wells tighter effectively makes the 
spin relaxation slower. 

This similarity between effects of increasing number of 
wells and increasing rso is not complete: one can see 
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FIG. 4. Dimensionless magneconductivity Sa(B) for 
Hso I H v = 2 and H\2 / H v — 4 and different number of wells 
TV = oo, 6, 4, 3, and N — 2 (solid lines, in order from higher 
to lower near the minimum). The solid lines were calculated 
without using the periodical boundary conditions. The mag- 
netoconductivity calculated with these conditions is shown for 
N = 4 by a dashed line. The inset shows the effect of spin re- 
laxation on magnetoconductivity for N = 2 and H12/ ' H v — 4. 

from Fig. |] that the minimum on the magnetoconduc- 
tivity curve does not disappear as N increases, while in- 
creasing rso leads to a monotonic 5a(B). After study- 
ing the dependencies of 5a(B) on all the parameters, we 
have found that if for a given rso and some N, T\i the 
function Sa(B) has a minimum, this will not change for 
any other N and Tn as long as rso is maintained the 
same. To better understand this behavior, we examine 
the magnetoconductivity dependence on t%2- Figure [5] 
shows the magnetoconductivity for N = 3, Hso /H v = 2, 
and Hi2/H v from to 100. From the analysis of the lim- 
iting cases T12 <C t v and T12 S> t v we can conclude that 

Sa(B, T12 = 00) = N5a(B, t X2 = 0). (30) 

It is obvious that both limiting cases are either mono- 
tonic, or not. From Fig. |^ one can see that, while it is 
not true that all 5a{B, 712) curves are contained between 
these two limits, the minimum is present on all of the 
curves. 

The situation at N — 00 is slightly different, since 
from the analysis of the limiting cases and Eq. ( ^0| ) it 
follows that the conductivity "per one well" 6a tends to 
at small t\i (the conductivity of the entire system is, 
of course, infinite for any T12 > 0). Therefore, with de- 
creasing T12 the magnetoconductivity also decreases, as 
demonstrated in Fig. |[ However, the main qualitative 
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FIG. 5. Dimensionless magnetoconductivity 5a for N = 3, 
Hso/H v = 2, and Hi2/H v from to 100. The squares show 
3 ■ Sa(B,Hi2/H ip = 100) (which is very close to the limit 
at H12 = 00), to illustrate the connection Eq. ( po| ) between 
magnetoconductivity at H12 = and H12 = 00. 

features of the magnetoconductivity dependence on tun- 
neling remain the same: if the minimum of 6 a exists 
with no tunneling, it also exists at any tunneling rate. 
The insert shows that the opposite is also true: if the de- 
pendence Sa(B) was monotonic for 712 = 00, it remains 
so for any finite t%2. 

In conclusion, we have analyzed the effect of weak 
tunneling between quantum wells in a multi-quantum- 
well structure on the weak localization and the mag- 
netoconductivity caused by it. We have found an an- 
alytical solution of the problem for the case when spin 
relaxation can be described by a single spin relaxation 
time, neglecting the effects of the linear terms in spin re- 
laxation Hamiltonian. The tunneling across the MQW 
structure was found to decrease the overall magnitude of 
the weak localization correction. However, it does not 
change the qualitative character of the magnetoconduc- 
tivity, namely, presence (or absence) of regions of posi- 
tive and negative magnetoconductivity. Our results show 
that the weak localization correction to the conductivity 
in multi-quantum-well structures is sensitive to the tun- 
neling time T12 and it should be possible to determine 
this time experimentally from the magnetoconductivity 
measurements . 

G. E. P. acknowledges support by RFFI Grant 96-02- 
17849 and by the Volkswagen Foundation. 
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FIG. 6. Dimensionless magnetoconductivity 8a for 
N = oo, H so /H v = 2, and H 12 /H v from to 100. The 
inset shows magnetoconductivity for Hso/H v = 0.5. 
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